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Abstract 

We present the microscopic equation for the growing interface with quenched 

noise for the model first presented by Buldyrev et al. [Phys. Rev. A 45, 

R8313 (1992)]. The evolution equation for the height, the mean height, and 

the roughness are reached in a simple way. The microscopic equation allows 

us to express these equations in two contributions: the contact and the local 

one. We compare this two contributions with the ones obtained for the Tang 

and Leschhorn model [Phys. Rev A 45, R8309 (1992)] by Braunstein et al. 

[Physica A 266, 308 (1999)]. Even when the microscopic mechanisms are quiet 

different in both model, the two contribution are qualitatively similar. An 

interesting result is that the diffusion contribution, in the Tang and Leschhorn 

model, and the contact one, in the Buldyrev model, leads to an increase of 

the roughness near the criticality. 
PACS numbers: 47.55.Mh, 68.35.Fx 
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I. INTRODUCTION 



The description of the noise-driven growth on a self-affine interface far from equilibrium is 
a challenging problem. The interface has been characterized through scaling of the interfacial 
width w with time t and lateral size L. The result is the determination of two exponents 
(3 and a called dynamical and roughness exponents, respectively. It is well known that 
interfacial width w ~ L a for t ^> t* and w ~ t@ for t <C t*, where t* ~ L a ^ is the saturation 
time. These properties occur for many models of surface growth. The values of the exponents 
leads to the classification of these models in different universality classes. Several models, 
belonging to the same directed percolation depinning (DPD) universality class, have been 
introduced to explain experiments on fluid imbibition in porous media, roughening in slow 
combustion of paper, growth of bacterial colonies, etc. 

It is currently accepted that the quenched disorder plays an essential role in those exper- 
iments. The DPD models take into account the most important features of the experiments 
|]J3|]. The two first models were simultaneously introduced by Buldyrev et al. [Q and Tang 
and Leschhorn || to explain the fluid imbibition in paper sheet. Several authors have been 
focused their attention on scaling properties and relationships between the dynamical and 
statics exponent for these models. The Tang and Leschhorn (TL) model has been recently 
reviewed by Braunstein et al. from a different point of view than the traditional one. 
The principal contribution was the restatement of the Microscopic Equation (ME) for the 
TL model (see Appendix). This equation allows the separation into two contributions: the 
substratum contribution by local growth and the diffusion one. They found that the diffu- 
sion contribution to the temporal derivative of the square roughness may be either negative 
or positive and that the behavior of this contribution depends on the pressure p. The neg- 
ative contribution tends to smooth out the surface, this case dominate for p > p c (where 
p c = 0.461 is the critical pressure). The positive contribution enhances the roughness. At 
the critical pressure the substratum contribution to the temporal derivative of the square 
roughness is practically constant, but the diffusion contribution is very strong. This last 
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contribution, has important duties on the power law behavior. 

In this paper we focus the attention in the Buldyrev et al. model of DPD. We show that 
this model presents several qualitative features of the TL model. We write a ME, starting 
from the microscopic rules, for the evolution of the height as a function of time. The ME 
allows us the identification of two contributions that dominates the dynamics of the system, 
the "contact" and the "local" one. In this context we study the mean height speed (MHS) 
and the temporal derivative of the square interface width (DSIW) as a function of time. We 
show that the contact contribution smooth out the surface for p well above the criticality 
but enhances the roughness near the critical value. To our knowledge the separation into 
two contributions for all the quantities studied in the present paper and the important 
duties of the contact contribution to the critical power law behavior has never been studied 
before. The paper is organized as follows. In Section |T| we write the microscopic equation 
for the evolution of height, the mean height and the roughness for the Buldyrev et al. 
model. We study the MHS, analyzing the contact and the local contributions. Also, the 
two contributions to the DSIW are analyzed. This separation into two contributions allows 
us to explain the mechanism of roughening. In Section |TJ we compare the Buldyrev et al. 
model with the TL model. Finally, we conclude with a discussion in Section fVj. 



II. BULDYREV ET AL. MODEL 

A. Microscopic Equation 

The interface growth takes place in a lattice of edge L with periodic boundary conditions. 
A random pinning force g(r) uniformly distributed in [0, 1] is assign to every cell of the 
lattice. For a given pressure p, the cells are divided in two groups, active (free) cells with 
g(v) < p and inactive (blocked) cells with g(r) > p. The interface between wet and dry cells 
is specified by a set of integer column heights hi (i = 1, . . . , L). At t = we start with flat 
initial conditions, i.e. hi = for all %. During the growth, a column is selected at random 
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with probability 1/L and the highest dry active cell, in the chosen column, that is nearest 
neighbor to a wet cell is wetted. Afterwards, we wet all the dry cell below it. In this model, 
the time unit is defined as one growth attempt. In numerical simulations at each growth 
attempt, the time t is increased by St = 1/L. In this way, after L growth attempts, the time 
is increased by one unit. In our simulations we use L = 8192. 

We consider the evolution for the height of the z-th site of the process described above. 
Let us denote by hi(t) the height of the i-th generic site at time t. From the microscopic 
rules we obtain the evolution for the i-th height in the next time step St = 1/L 

tuft + st) = Kit) + st {e{- Zi ) + i) + [i - e(-zi)) Yi} , (i) 

where Q(x) is the unit step function defined as Q(x) = 1 for x > and equals to 
otherwise, Fiihi + j) equal to 1 if the cell is active and if the cell is inactive (1 < j '• < Zi), 
Zi = max(/ij_i, hi+i) — hi, and 

Yi = ZiFi(hi + z^ 

+ ( Zi - 1) Fi(hi + Zi - 1)(1 - Fi(hi + z^) 

+ {zi - 2) Fi(hi + Zi-2)(1- Fi(hi + Zi - 1)(1 - Fi(hi + Zi)) + ■■■ 
+ Fi(hi + 1)(1 - Fi(hi + 2)) ... (1 - Fi(hi + Zi - 1)(1 - Fi(hi + Zi)) , 

is the increase of the height in the i-th column due to the contribution of the nearest lateral 
neighbor. The term between braces in Eq. (|l]) takes into account all the possible ways the 
site % can growth. The height in the site % is increased by 

1. 1 if hi > max(/ij +1 , h^i) and i^(/tj + 1) = 1, 

2. Yi if hi < max(/tj + i, hi-i) and + Y/) = 1. 

Otherwise, the height is not increased. We shall call contact contribution to the term 
[1 — 0(— Zi)] Yi (related to case 2) and local contribution to the term Q(—Zi) Fi(hi + 1) 
(related to case 1). 
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Averaging Eq. ([!]) over the lattice, taking St — > the evolution equation for the mean 
height is 

^ = (e(-^)F < > + <(i-e(-^))y i >, (2) 

and the evolution equation for the square interface width is 
dw 2 



2({h i - (hi))e(-Zi) Fi) + 2{(h i - (hi))(l - e(-Zi)) Yj) . (3) 



dt 

The first terms of both equations can be identified as the local growth contribution, and the 



second term as the contact growth contribution. We shall see in Subsection [II Q that the 
separation into these two analytical terms allows us to show how the contact mechanism 
enhances the roughness near the criticality. In the present paper we focus only on the 
dynamical behavior, i.e. t <C t* ~ L for the mean height and roughness (in the DPD models 
a/P= 1). 



B. Mean Height 

The top plot of Figure |l| shows the MHS as a function of time in three regimes (moving, 
critical and pinning phases). The initial condition for the MHS is p in all regimes. At the 
criticality we found for the mean height a power law behavior with the same dynamical 
exponent that the roughness one (3 = 0.68 ± 0.02 for p c = 0.531. In the moving and 
pinning phase we can see that this power law does not hold. Bellow the criticality, in the 
pinning phase, the MHS go to zero. Above, in the moving phase, the MHS goes to certain 
constant value. The left plots of Figure show the contributions to the MHS: the local 
one (Q(-Zi)Fi) and the contact one ((1 — 0(— zi))Yi). The local contribution, which is 
equal to p at t — 0, is stronger in the early time regime. This is because in this regime the 
difference of heights between nearest neighbors is mostly less than one. This contribution 
set into motion the contact growth. In the moving phase (see left top plot of Figure 0) both 
contributions go to a certain constant. In the intermediate regime the local contribution 
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decreases while the contact one increases. At the criticality and in the pinning phase (see 
left middle and left bottom plots of Figure |2], respectively) the local contribution decreases 
continuously from p to zero. The contact contribution increases from zero to a maximum 
value and then decreases reaching asymptotically the MHS. In all phases both contributions 
are equal at t ~ 1. This means that after L-growth attempts the interplay between both 
mechanisms are equal independently of p. After this, the dynamical behavior is strongly 
dominated by the contact mechanism. 

C. Roughness 

The top plot of Figure |3] shows the temporal (DSIW) as a function of time for various 
values of p. The initial condition is p in all regimes. As we expected 0, the power law holds 
only at the criticality. The DSIW goes asymptotically to zero at the pinning and moving 
phase. In the left plots of Figure |], we show the two contributions to the DSIW for different 
values of p. The local contribution 2((hi — (/ij))6(— Zj) Fi) to the DSIW is always positive. 
As p decreases, this contribution becomes less important, but always rough the interface. 
On the other hand, for p > p c , the contact contribution 2((/ij — — G(— Zi)) Yj) can take 

negative values, smoothing out the surface. Otherwise, for p < p c , the contact contribution 
is always positive roughening the interface. One could expect that the contact contribution 
always smooth out the surface because it tends to widen the roughen picks. However, near 
the criticality, the contact growth happens mainly in lateral neighbors cells to few height 
terraces above the mean height. Then, this new wetted column smooth out locally, but it 
moves away from the mean height increasing the roughness. 

III. COMPARISONS WITH THE TANG AND LESCHHORN MODEL 

We rescue the similarities between the Buldyrev et al. and the Tang and Leschhorn 
models despite the strong microscopic differences between their rules. In a previous paper 



6 



Braunstein et al. @|| wrote the ME for the TL model. They identified two separate contri- 
butions: the substratum and the diffusion one in the MHS and the DSIW (see Appendix). In 
the present paper, we also obtain two contributions: the local and the contact one. Figure |2| 
shows the contributions to the MHS for both models. Notice that each pair of plots have the 
same qualitatively behavior. The shape of the diffusion and substratum contribution in the 
TL model are qualitatively similar with the contact and local contribution in the Buldyrev 
et al. model, respectively, even when the microscopic processes are quite different for each 
model. The different contributions to the DSIW for both models are shown in Figure f|. No- 
tice that the diffusion and the contact contributions, in each model, can take negative values 
for p > p c smoothing out the interface. Near the criticality, in both models, the roughness 
is mainly due to the diffusion and the contact contributions. These last contributions play 
a very important role at the criticality in each model. The similarities between the models 
could explain why these two different microscopic models belongs to the same universality 
class. Figure [3] shows the DSIW as a function of time for both models. As we expected ||, 
the power law holds only at the criticality despite other authors [0,0,|7J-|^] . From these plots 
it is easy to see that this last statement holds for both models. 



IV. CONCLUSIONS 

We wrote the ME for the evolution of the height in the Buldyrev et al. model and we 
compared the results obtained with those from the TL model. Using the ME we studied the 
evolution of the mean height and the roughness. The ME allows us to separate the local 
and the contact contributions. We found that the contact contribution near the criticality 
is the main responsibility of the roughness. We found qualitatively that the shape of the 
contact contribution is analogous to that of the diffusion one in the TL model, and that the 
shape of the local growth is similar to that of the substratum contribution in the TL model. 
We found that the power law behavior holds only at the criticality for the Buldyrev et al. 
model. This last feature, common to the TL model, suggests to us that it may be common 
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to all other DPD growth models 0. 
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APPENDIX: MICROSCOPIC EQUATION FOR TANG AND LESCHHORN 

MODEL 



We present here the microscopic equation for the TL model @J§]. The time evolution 
equation for the interface height, in a time step St, is 



tn(t + St) = h^t) + St [W i+1 + Wi-i + F^K) Wi] , (Al) 



with 



Wi = l- Q{hi - min(/ii_i, h i+l ) - 2) , 
Wi±i = Q(h i±1 - mm(hi, h i±2 ) - 2){[1 - 0{hi - h i±2 )} + \S hiM±2 } . 

where — hi + 1. Here Q(x) is the unit step function defined as Q(x) = 1 for x > and 
equals to otherwise. Fi{h'j) equals to 1 if the cell at the height h\ is free or active (i.e. the 
growth may occur at the next step) or if the cell is blocked or inactive. F^ is the interface 
activity function. 

Averaging over the lattice, taking St — > the evolution equation for the mean height is 

^ = (W i F i ) + {l-W i ), (A2) 

and the evolution equation for the square interface width is 

— = 2((hi - (h^WiFi) + 2<[min(/i i _ 1 , h i+1 ) - (h^l - W { )) . (A3) 

The first terms of both equations have been identified as the substratum growth contribu- 
tions, and the second terms as the diffusion growth contribution. 
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FIGURES 

FIG. 1. Plots oip~ l dh/dt vs t. The top (bottom) plot shows the results for the Buldyrev et 
al. (Tang and Leschhorn) model. For the top plot the parameter p is 0.56 (A), 0.531 (O) and 0.51 
(y). For the bottom plot the parameter p is 0.49 (A), 0.461 (O) and 0.4 (y). 

FIG. 2. ln-ln plots of different contributions to the MHS as a function of time for the Buldyrev 
et al. model (left plots) and the Tang and Leschhorn model (right plots), for different values of 
p. The circles represent the contact contribution (left plots) and the diffusion contribution (right 
plots). The squares represent the local contribution (left plots) and the substratum contribution 
(right plots). For both models, the top, bottom and middle plots show the moving, the pinning 
and the critical phases, respectively. 

FIG. 3. DSIW as a function of time. The top (bottom) plot shows the results for the Buldyrev 
et al. (Tang and Leschhorn) model. For the top plot the parameter p is 0.7 (0)> 0.56 (y), 0.531 
(•) and 0.51 (A). For the bottom plot the parameter p is 0.7 (O), 0.49 (y), 0.461 (•) and 0.4 
(A). In both models the symbol • shows the critical behavior. 

FIG. 4. Semi- In plots of different contributions to the DSIW as function of time for the Buldyrev 
et al. model (left plots) and the Tang and Leschhorn model (right plots), for different values of 
p. The circles represent the contact contribution (left plots) and the diffusion contribution (right 
plots). The squares represent the local contribution (left plots) and the substratum contribution 
(right plots). For both models, the top, bottom and middle plots show the moving, the pinning 
and the critical regime, respectively. 
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